THE QUANTUM HYDRODYNAMICS SYSTEM IN TWO SPACE 

DIMENSIONS 



g . PAOLO ANTONELLI AND PIERANGELO MARCATI 

' Abstract. In this paper we study global existence of weak solutions for the 

^ ' Quantum Hydrodynamics System in 2-D in the space of energy. We do not 

require any additional regularity and/or smallness assumptions on the initial 
data. Our approach replaces the WKB formalism with a polar decomposition 
f**^ theory which is not limited by the presence of vacuum regions. In this way 

' we set up a self consistent theory, based only on particle density and current 

density, which does not need to define velocity fields in the nodal regions. 
The mathematical techniques we use in this paper are based on uniform (with 
' respect to the approximating parameter) Strichartz estimates and the local 

' smoothing property. 

1. Introduction 

In this paper we are concerned with the Cauchy problem for the Quantum Hy- 
drodynamics System in two space dimensions 

^ r %>+divj = o 

(i.i) I a t j + div(^)+vp(p) + P v^ + /(vp,J,vvp) = fpv(^) 



in 



-AV = p-C(x), 
where t > 0, x £ M 3 , with initial data 
(1-2) p(0)=po, J(0) = Jo- 

The unknowns p, J represent the charge and the current densities, respectively, 
P(p) is the classical pressure, which we assume to satisfy the usual 7— law, P(p) = 
^pj-pT - , with p to be specified later. V is the self-consistent electric potential, 

given by the Poisson equation. The term ^-pV y ^fp j can t> e interpreted as the 

quantum Bohm potential, or as a quantum correction to the pressure. Indeed with 

some regularity assumption we can write the dispersive term in different ways: 

fr 2 /A fo\ fr 2 fr 2 
(1.3) y pV(^-^J = T div(pV 2 lo g(0 ) = T AV P -^div(VVp®VVp). 

There is a formal analogy between (1.1) and the classical fluid mechanics system, 
in particular when h = 0, the system (1.1) formally coincides with the (nonhomo- 
geneous) Euler-Poisson incompressible fluid system. 

The theoretical description of microphysical systems is generally based on the wave 
mechanics of Schrodinger , the matrix mechanics of Heisenberg or the path-integral 
mechanics of Feynman (see [12]). Another approach to Quantum mechanics was 
taken by Madelung and de Broglie (see [23]), in particular the hydrodynamic theory 
of quantum mechanics has been later extended by de Broglie (the idea of "double 
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solution") and used as a scheme for quasicausal interpretation of microphysical 
systems. There is an extensive literature (see for example [10], [21], [19], [18], and 
references therein) where superfluidity phenomena are described by means of quan- 
tum hydrodynamic systems. Furthermore, the Quantum Hydrodynamics system 
is well known in literature since it has been used in modeling semiconductor de- 
vices at nanometric scales (see [13]). The hydrodynamical formulation for quantum 
mechanics is quite useful with respect to other descriptions for semiconductor de- 
vices, such as those based on Wigner-Poisson or Schrodinger-Poisson, since kinetic 
or Schrodinger equations are computationally very expensive. For a derivation of 
the QHD system we refer to [1], [17], [14], [11]. 

In (1.1) the term /(^/p, J, V^/p) represents collisions between electrons in models 
for semiconductor devices, see for example [4] where f = J. Otherwise, more gen- 
erally it is used to mimic the interaction of the quantum fluid with an external 
fluid. In this paper we will treat the case / = J, whereas we will consider the more 
general case in a forthcoming paper. 

We are interested in studying the global existence of weak solutions to the system 
in the class of finite energy initial data without higher regularity hypotheses or 
smallness assumptions. Local existence of solutions to (1.1) is proved in [16], while 
global existence with further regularity for the initial data is proved in [22]. For 
the three dimensional case the problem of proving existence of weak solutions in 
the space of energy is solved in [2]. 

There is a formal equivalence between the system (1.1) and the following nonlinear 
Schrodinger-Poisson system 



(1.4) 



ihd t ip + ^-Aip = |V| P_1 ^ + Vip + Vi) 
-AV = \ip\ 2 

where V = ^ log (jt^J , in particular the hydrodynamic system (1.1) can be obtained 

by defining p — \ip\ 2 , J = him (ipVip) and by computing the related balance laws. 
This problem has to face a serious mathematical difficulty connected with the need 
to solve (1.4) with the ill-posed potential V. 

On the other hand, in two space dimensions we have to face another mathematical 
difficulty, which is the singularity of the electrostatic potential V. Indeed, the 
Green's function for the Poisson equation in R 2 reads 

(1.5) *(ar) := -JLlog|a:|, x e K 2 , 

and consequently the electrostatic potential is implicitly given by the following 
convolution 

(1.6) V(t,x) := --!- / log | a; - y\p(t,y)dy. 

Hence V could be unbounded at infinity, even if we assume p to be sufficiently 
decaying at infinity. Consequently some additional hypotheses on the electrostatic 
potential at the initial time are required in order to assure its boundcdness and 
some suitable integrability properties. 

A natural framework to study the existence of the weak solutions to (1.1) is given 
by the space of finite energy states. Here the energy associated to (1.1) is defined 

by 

(1.7) E(t) := y |V7^)| 2 + i|A(t)| 2 + f(p{t)) + i|W(t)| 2 dx, 

where A := J/y/p, and f(p) = ^jP^ 1 '■ The function f(p) denotes the internal 
energy, which is related to the pressure through the identity P{p) = pf'(p) — f(p)- 
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Therefore our initial data are required to satisfy E(0) < oo and to have finite mass, 
i.e. / podx < oo, or equivalently (by the Gagliardo-Nirenberg inequality), they are 
required to be such that 

(1.8) V^eff^R 2 ) and A„ := J /^e L 2 (t 2 ). 

This paper is organized as follows: in Section 2 we present the main theorem, 
and the notations and main mathematical tools used throughout the paper. In 
Section 3 we recall the polar decomposition of a wave function into its amplitude 
and a unitary factor, as presented in [2]. In Section 4 we give a result of global 
existence of solutions to a nonlinear Schrodinger-Poisson system in two dimensions, 
and explain its equivalence with the QHD system (1.1) without a collision term. In 
Section 5 we present how to construct a sequence of approximate solutions to (1.1), 
and show the consistency for this approximate solutions. In Section 6 we prove 
some a priori bounds for the sequence of approximate solutions in order to pass 
to the limit, prove the convergence and showing the existence of weak solutions to 
(1.1). 

2. Statement of the main theorem and preliminaries 

First of all, let us recall the definition of a weak solution to (1.1), (1-2). 

Definition 2.1. We say the pair (p, J) is a weak solution to the Cauchy problem 
(1.1), (1.2) in [0,T) x R 2 with locally integrable initial data (po, Jo), if there exist 
VP, A, such that VP G L 2 oc ([0,T) ;j ffi oc (R 2 )), A e L 2 OC ([0, T); L 2 oc (R 2 )) and by 
defining p := (yfp) 2 , J := ^/pA, one has 

• for any test function r\ e Cg°([0,T) x R 2 ) we have 

(2.1) / / pd t rj + J ■ \7rjdxdt + f p r](0)dx = 0; 
Jo Jr 2 Jr 2 

• for any test function C € Cg°([0,T) x R 2 ;R 2 ) 

(2.2) / / J-d t ( + A®A: V( + P(p)div(- pVV ■(- J ■ ( 
Jo Jr 2 

h 2 f 
+ h 2 Vi/p <g> VVp" : VC - -rpA div Qdxdt + Jo ■ ({0)dx = 0: 

4 J R 2 

• (generalized irrotationality condition) the identity 

(2.3) VAJ = 2VVpAA 
holds in the sense of distributions. 

Remark 2.2. Suppose we are in the smooth case, so that we can factorize J = 
pu, for some current velocity field u, then the last condition (2.3) simply means 
pV Am = 0, i.e. the current velocity u is irrotational in pdx. This is why we will 
call it generalized irrotationality condition. 

Definition 2.3. We say that the weak solution (p, J) to the Cauchy problem (1.1), 
(1.2) is a finite energy weak solution (FEWS) in [0, T) x R 2 , if in addition the energy 
(1.7) is finite for almost every t € [0,T). 

The hydrodynamic structure of the system (1.1), (1.2) should not lead to con- 
clude that its solutions behave like classical fluids, indeed the connection with 
Schrodinger equations suggests that (p ,J ) should in any case be seen as mo- 
menta related to some wave function tpQ. The main result of this paper is the 
global existence of FEWS by assuming only the initial data (p , Jo) are momenta 
of a wave function ip e H 1 (R 2 ), without smallness or smoothness additional con- 
ditions. Furthermore, as already pointed out above, because of the singularity of 
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the electrostatic potential, we have: to make an additional assumption on the initial 
mass density po. We can now state the main theorem of this paper. 

Theorem 2.4. For any ip £ H 1 (R 2 ) let us denote 

Po = IV'oI 2 , = hIm(ip Vip ), 

and furthermore assume that 

(2.4) / p a \ogp dx < oo, 
and 

(2.5) V {0,x):=~ f log |a; - y\p (y)dy £ 17 (E 2 ), 

2tt J R 2 

for some 2 < r < oo. Then there exists a FEWS of (1.1), (1.2), in [0,T) x K 2 . 

2.1. Notations. For convenience of the reader we will recall some notations which 
will be used in the sequel. 

If X, Y are two quantities (typically non- negative), we use X < Y to denote 
X < CY, for some absolute constant C > 0. 

We will use the standard Lebesgue norms for complex- valued measurable functions 
/ : R d C 



\Lp( 



\f(x)\Pdx 



If we replace C by a Banach space X, we will adopt the notation 

Wfh^x) ■= (Jjimwxdx^) /P 

to denote the norm of / : R d — > X. In particular, if X is a Lebesgue space L r (R n ), 
and d = 1, we will shorten the notation by writing 



< K n) := ll/Wllir^jdt) ^ - (jT^ |/(t ;a; )rdx)^di 



1/q / r r x 1/9 

L«Lj(/xI 

to denote the mixed Lebesgue norm of / : I — > L r (M. n ); moreover, we will write 
L\L r x {I x E") := Li(I; L r (E™)). 

For s e E wc will define the Sobolev space H s (W l ) := (1 - A)- S / 2 L 2 (M"). 
2.2. Nonlinear Schrodinger equations setting. 

Definition 2.5 (see [8]). We say that (q,r) is an admissible pair of exponents if 
2<g<oo,2<r<oo, and 

, . 1 1 1 

^ - q + r = 2- 

Now we will introduce the Strichartz norms. For more details, we refer the reader 
to [26]. Let I x K 2 be a space-time slab, we define the Strichartz norm S°(I x M 2 ) 

IMIs°( 7 X R2 ) : = SU P ll u IU?L S (/xR2), 

where the sup is taken over all the admissible pairs (q,r). For any fc > 1, we can 
define 

ll u lls ,fc (/xR 2 ) := \\^ ku \\s°(IxM 2 )- 
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Theorem 2.6 (Keel, Tao [20]). Let (q,r), (q,f) be two arbitrary admissible pairs 
of exponents, and let U(-) be the free Schrddinger group. Then we have 

(2-7) \Mt)f\\ L , L r < \\f\\ L 2 m 

(2-8) || / U(t-s)F( S )ds\\ L <, L r<\\F\\ Li , L „ 

J s<t ' x 

(2-9) || J U(t-s)F(s)ds\\ L , L r<\\F\\ LfLf ,. 

By using the previous Theorem, we can state the following Lemma about some 
further integrability properties for solutions of the Schrodinger equations. 

Lemma 2.7. Let I be a compact interval, and let u : I x R 2 — > C be a Schwartz 
solution to the Schrddinger equation 

idtu + Au = F 1 + .. . + F M , 

for some Schwartz functions F\ , . . . , Fm ■ Then we have 

IMWxr") < \\u(to)\\ w + ll^ll^^ (7xR2) + • • • + II^ii l: ^ (/xR2) , 

where (qi, ri), . . . , (qM, I'm) are arbitrary admissible pairs of exponents. 

Furthermore, the solutions to the Schrodinger equation enjoy some nice local 
smoothing properties; there are various results in the literature regarding this prop- 
erty, here we recall a theorem, due to Constatin and Saut [9] which actually cov- 
ers a more general setting. In particular it can been shown the solutions of the 
Schrodinger equation are locally more regular than their initial data. Close results 
have been proved by P. Sjolin [25] and L. Vega [29]. 



Theorem 2.8 (Constantin, Saut [9]). Let u solves the free Schrddinger equation 

id t u + Au = 
w(0) = uq. 



(2-10) 



Let X € Of (R 1+2 ) of the form 

X(t,x) = Xo(t)xi(xi)X2(x2), 
with Xj € Co° W • Then we have 

f X 2 {t,x)\{L - Ay/ 4 u(t,x)\ 2 dxdt < C 2 \\u Q \\ LH 

In particular, if u € L 2 (R 2 ), one has for all T > 

U a 2 ([0,T];<(I 2 )). 
We have a similar result also for the nonhomogeneous case: 
Theorem 2.9 (Constantin, Saut [9]). Let u be the solution of 

(2.H) 

where F G L 1 ([0, T]; L 2 (R 2 )). Then it follows 

ueL 2 ([0,T];HH 2 (R 2 )). 
Moreover, let x G Cg°(M 1+2 ) be of the form 

X(t,x) = Xo(t)xi(xi)X2(x 2 ), 



id t u + Au = F 
u(0)=u o G L 2 (R 2 ), 
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with Xj € suppxo C [0,T]. Then the following local smoothing estimate 

holds 

(2.12) (f X 2 {t 1 x)\{I-K) 1 / i u{t,x)\ 2 dxdt\ 

< C (JI u o||l 2 (r 2 ) + II^IIlJl^io.tIxr 2 ) ■ 

This results imply that, the free Schrodingcr group U(-) fulfills the following 
inequalities 



(2.13) ||tf(>o|| 

(2-14) || J U(t- s)F(s)ds\\ L2{[0 T] . H i/2 {R3)) < \\F\\ L i mT] . L 2 (m} . 

2.3. Logarithmic Sobolev inequality and the Gagliardo-Nirenberg inequal- 
ity in 2D. As we already pointed out in the introduction, in two dimensions the 
Poisson electrostatic potential is highly singular, because of its kernel. Hence the 
need to make the further assumption on the initial datum p , as stated in Theorem 

2.4. A fundamental tool to ensure for the energy to be bounded at all times is the 
logarithmic Sobolev inequality. 

Theorem 2.10 (Carlen, Loss [7]). Let f be a non-negative function in L X (R 2 ) such 
that /log/,/log(l + |.t| 2 ) e L^R 2 ). Ifjfdx = M, then 

Mr r 

(2.15) — / /log /da: + / /(*) log \x - y\f(y)dxdy > 

1 JR 2 JR 2 xR 2 

M 2 

> C(M) := — (l + log7r + logM). 

Furthermore in the sequel we use the classical Gagliardo-Nirenberg inequality in 
R 2 , namely 

(2.16) hiiL P <hii! 2 iiv U |ii; f , 

for 2 < p < oo. 

2.4. Compactness tools. In conclusion we recall a compactness result due to 
Rakotoson, Tcmam [24] in the spirit of classical results of Aubin, Lions and Simon 
which will be needed in Section 6. 

Theorem 2.11 (Rakotoson, Temam [24]). Let (V, \\-\\v), (H; \\-\\h) be two separable 
Hilbert spaces. Assume that V C H with a compact and dense embedding. Consider 
a sequence {u e }, converging weakly to a function u in L 2 ([0,T]; V), T < oo. Then 
u e converges strongly to u in £ 2 ([0, T]; H), if and only if 

(1) u £ (t) converges to u(t) weakly in H for a.e. t; 

(2) lim| EKO ,.Ec[o ) T] su Pe>0 J E Hu^H^dt = 0. 

3. Polar decomposition 

In this Section we will recall (see [2]) the main properties related to the polar 
factorization of wave functions (detalis will be given in the Appendix). The polar 
decomposition is a useful tool to factorize a wave function ip into its amplitude y/p = 
\ip\ and a unitary factor <f> such that ip = y/p4>, that we write the hydrodynamical 
quantities V-s/p, A without the necessity of defining the velocity u = J/p and the 
phase S as in the WKB formalism. 

Let us consider a wave function ip <G L 2 (R 2 ) and define the set 

(3.1) P(ip) '■= {4 1 measurable | ||<^||loo(r2) < 1, yfpcf) = ip a.e. in R 2 }, 
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where yfp = \xp\. Of course if we consider <f> £ -P(V0> then by the definition of 
the set P{ip) it is immediate that |</>| = 1 a.e. in ^fpdx and <fi G P(ip) is uniquely 
determined a.e. ^fpdx in R 2 . 

Through the polar factorization of the wave function ip we get the following identity 
holds 

7 63 T 

(3.2) ft 2 Rc(W <8> WO = — y- + fi 2 VVP ® V yfp. 

Indeed by writing the hydrodynamical quantities V-s/p, A in terms of the associated 
wave function tp and its polar factor <fi, it follows 

(3.3) VVp = Rc(0VV>), A = Mm(^W). 

Hence the equation (3.2) holds true since it hides a null structure, namely 

h 2 Re{((j)VW) ® (0Vr/O} = ?X 2 Re(^V^) <g> Re(0VV>) + h 2 Im(0Vip) <g> Im(^W). 

Moreover this bilinear structure is H 1 — stable, namely if {tp n } C ff 1 (M 2 ) is a 
strongly convergent sequence in H , ^> n — > -0, then also their related hydrody- 
namical quantities converge strongly, Vy/pn — > V^//3, A„ — > A in i? 1 (R 2 ). This 
is a quite remarkable property, since for the polar factors we could only have a 
weak-* convergence in L°° (up to passing to a subsequence), since as we already 
pointed out they are uniquely defined only in ^fpdx. Hence the weak convergence 
of the quantities Re(<fi n Vip n ), lm((f) n Vip n ) is upgraded to the strong one, thanks to 
the bilinear structure (3.2). Furthermore, it is worth noting that this structure is 
fundamental for the stability in the energy space: with the WKB formalism for 
example this couldn't be achieved without some further regularity assumptions. 

Lemma 3.1 (Existence and Stability Lemma). Let tp <G H 1 (M. 2 ), y/p := \ip\, then: 

(i) there exists <j> g P(tp) C L°°(IR 2 ) such that tp = Jpfy a.e. inR 2 , ^fp g 
H 1 (R 2 ), (f> is unique in ^fpdx, Vy 7 /? = Re(^VV'), if we denote A := 
Mm($Vip), then A g L 2 (R 2 ) and moreover 

(3.4) h 2 Re(d^d k iP) = h 2 d jy /pd ky /p + A (j) A w 
and 



(3.5) hVip A V?/> = 2iVy/p A A. 

(ii) (Stability) Let ip n — > ip strongly in H 1 (R 2 ), denote by <j) n the polar factor 
ofipn, then it follows 

(3.6) VVP^VVP, A„^A mL 2 (R 2 ), 
where A„ := hIm(<f) n S7ip n ). 

The previous identity (3.5) is needed to obtain the generalized irrotationality 
condition (2.3), since V A J = Mm(Vi/> A Vip). 

In the Bohmian mechanics (see [28]) and in the Nelson's stochastic mechanics (see 
[6]) the major difficulty regards the definition of current and osmotic velocities in 
the nodal regions (p = 0). The main advantage of our approach is to build up a 
self consistent theory, based only on particle density and current density, which, 
thanks to the polar decomposition, does not require the use of velocity fields. 
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4. QHD WITHOUT COLLISIONS AND NONLINEAR SCHRODINGER-POISSON 

SYSTEM IN 2D 

In this section we consider the collisionless case, namely / = 0, 

!d t p + div J = 
d t J + div (^) + VP(p) + pVV = ^p\7 (^f) 
-AV = p, 

and our goal is to show the connection with the nonlinear Schrodinger-Poisson 
system 

{ } \-AV = |Vf, 

with initial datum "0(0) = ipo. 

Indeed, let us consider a (global) solution to (4.2), and the equation for the mass 
density p — \tp\ 2 is immediately given by 

d t p + div J = 0, 

where J = hlm^Vip). Now, if we compute the time derivative of the current 
density J, we find out the following conservation law 

ti 2 

d t J + div (h 2 Re(Vip <8> W>)) - — AVp + VP{p) + pVV = 0, 

where V is the Poisson electrostatic potential given by — AV^ = p, and P{p) = 
£_-Ip(p+ 1 )/ 2 . Now, from formula (3.4) wc can write 

h 2 Re(V^p ® W) = h 2 V^P ® VVp + A <g> A, 

where A is defined as in Lemma 3.1. Finally, we just notice that we can write, as 
in equation (1.3) 

^AVp - div(VVp® vvp) - \ P v {^jf) . 

where the distribution T = i/oV (^r) is defined by 

(T,<t>)=[ Vy/p^Vy/p :V(t>- ]-p A div cpdx, forall^eCS°(R 2 ), 
so that now the conservation law for the current density read 

(4.3) d t J + div (^) + VP(p) + pVF = f pV (^) . 

Hence, given a solution ^ to (4.2), with initial datum V(0) = ipo, by defining 
(p, J) := (\ip\ 2 , hlm(ipWip)) we have a solution to (4.1) with initial data (po,Jo) = 
(IV'ol 2 ! ^-I m (V'oVV'o))- These arguments can be written in a rigorous way by using 
a regularization argument as in Proposition 15 in [2]. 

Proposition 4.1. Letip <= H 1 (R 2 ), define the initial data, (po,Jo) '■= (l^ol 2 , ftlm(-0 o V^o)); 
moreover assume the following conditions hold: 

(4.4) / p \ogp dx < oo, 

Jr 2 

and 

(4.5) y(o,.):=-J- /" log|.-y|p (y)dyeL r (R 2 ), 

2tt J R 2 
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for some 2 < r < oo. Then for any T > there exists a FEWS (p,J) to the 
Cauchy problem (4.1) in the space-time slab [0, T) x R 2 . Furthermore the energy 
E(f) defined as in (1.7) is conserved for all times t £ [0,T). 

Consequently, to prove the Proposition 4.1 we only need to state a global well- 
posedness result for the Cauchy problem associated to equation (4.2). Indeed, let 
us point out that conditions (4.4) and (4.5) are needed to provide a global solution 
to the Cauchy problem associated to equation (4.2). Furthermore, as we will see in 
the next section, we also need the solution to (4.2) to satisfy appropriate estimates 
in some space-time mixed Lebesgue spaces. The result we prove is the following 

Theorem 4.2. Let us consider the Cauchy problem (4.2) with initial datum 

(4.6) V(0) = V>o, 
such that 

(4.7) Voe-H" 1 ^ 2 ), / |V-o| 2 log|V>o| 2 dx<oo, 

Jr 2 

and 

(4.8) y(0,.) : =--L f \og\--y\p (y)dy£L r (R 2 ), 

for some 2 < r < oo. Then, there exists a unique solution ip £ C([0, oo); iJ 1 (IR 2 )) 
to (4.2), (4.6). Furthermore the energy 

(4.9) E{t):=f ||VVf + ^h/T +1 + ^| 2 d* 

J R 2 2 p+1 2 



/ 



'j-\vyj( t ,x)\ 2 + ^m,x)\p +i dx 



p(t, x) log \x - y\p(t, y)dxdy, 



Ip 2 p + 

1 

47T 

is well-defined for all times and it is conserved. Finally, for every T > and for 
every (q,r) admissible pair of Schrddinger exponents in R 2 , we have 

(4-10) ||VV>||l?lj([o,t)xr*) < C(T, ||Vo||l',£o). 

Proof. Let us remark that now we're going to deal with the case in which the power- 
like nonlinearity |t/j| p-1 ?/> satisfies 1 < p < oo; the case p = 1 is straightforward. 
Firstly, from the logarithimic Sobolev inequality (see theorem 2.10), we infer the 
boundedness of the initial energy, and furthermore, since it is conserved for all 
times, we have 

E{t)=E < \\Mh + Wo\\h+ J Pologpoda;<oo. 

Secondly, since we have ip £ H 1 (M 2 ), from the Gagliardo-Nirenberg-Sobolev in- 
equality (2.16) we can say that H^Hl^r 2 ) for each 2 < p < oo, and consequently 
p £ LP(M. 2 ) for each 1 < p < oo. 

Now, we want to give a local well-posedness result for the solutions to (4.2). From 
standard arguments (see [26]), we look for a priori estimates of solutions to (4.2) 
in the Strichartz spaces. To this aim, we first need to know the mixed space-time 
Lebesgue spaces in which the electrostatic potential lies. 

Lemma 4.3. Let (ip, V) satisfy the S chro ding er-Pois son system (4.2) in [0, T) x R 2 
for some T > 0. Let us further assume that V(0) £ L r (R 2 ), for some 2 < p < oo. 
Then 

f WGL-LP([0,T)xR 2 ), Vpe[2,oo), 
1 ' \ V £ LfL r x ([Q,T) x R 2 ). 
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Remark 4.4. Let us note that the first a priori estimates in (4.11) for the gradient 
of the electrostatic potential don't imply any estimates on the potential V, since 
we can't use the Sobolev embedding theorem. Hence we need some additional 
assumptions on V^(O), namely hypothesis (4.8). 

Proof. Clearly by the elliptic regularity theory and by the Sobolev embedding the- 
orem we have 

(4.12) nwiUoo ia < mvrvii™ < lipoid- 

Furthermore, by taking the gradient of (1.6) we get 

(4.13) ^ V {t,x) = ~f j^Lp(t,y)dy, 

2tt J R 2 \x - y\ z 

and by using the generalized Hardy-Littlewood-Sobolev's inequality we obtain that 
W e L°°([0, T); L P (R 2 )), 2 < p < oo. Resuming, we infer W G L°°([0, T); L P (M 2 )), 
for 2 < p < oo. 

Let us now differentiate with respect to time the formula (1.6), we thus obtain, by 
using the continuity equation, 

dtV{t,x) = ~f . J(t,y)dy, 

2ir J R2 \x - y\ 2 

and consequently 

(4.14) V(t, x) = V(0, x) - 1- f f . j( S: y ) dyds . 

2tt Jo Jr 2 F - y\ 

Now, from the bounds on the energy and the Gagliardo-Nirenberg-Sobolev inequal- 
ity we infer J e L°°([0, oo); L P (R 2 )), for 1 < p < 2, thus the generalized Hardy- 
Littlewood-Sobolev inequality yields d t V e L°°([0, oo); L r (R 2 )), with 2 < r < oo. 
Hence 

(4.15) |in*)IU'<ll^(0)||Lr+t||ft^|| L?i:([0iOo)xR2)) 2<r<oo, 

and consequently V G L^°L r ([0, T) x M 2 ). □ 

Let us now come back now to the proof of theorem 4.2. By Strichartz estimates 
we have 



L t 1 L x 1 ([0,T)xR 2 ) 



+ ll^ll^^ ([ 0,T)x R2 ) + ll W ^ll^^ ([0>T)xR2)) 

where (q^r^) are dual exponents of an admissible pair for Schrodinger in R 2 , i 
1,2,3. Let us consider the second term in the right hand side, |||V'I P_1 W>|| „/ r 

By Holder's inequality we have 



where we have the following algebraic conditions on the exponents r[ , r 



< i < 



i<'4 = i + i<i 

2 — r r T 2 ^ 1 



Then, by applying the Gagliardo-Nirenberg (2.16) inequality we get 
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For the term || VV^>|| q ' r ' 2 we choose (g^r^) = (r 1 , 7x5), where r is the exponent 
such that V(0) € L r , and hence 

Hw^ii ^ <^||F|| L ~ iE ||vviU r ^ 

<T^(||V(0)||Lr+T||ftV||Loo/Ir)||V^|Ucc iS . 

Finally for the last term || VV^|| ,j r j we can easily choose three exponents r' 3 ,p,pi, 
such that 



and 

U ^ p - 2> U ^ pi - 2 

±<^ = ± + ^<l. 

2 — P ' Pi 

Hence we obtain the following estimate 

(4.16) ||VV||i«£r ( [o,DxR*) £ HVVolU- +r Q C(|^o||^,^), 

for some a > 0. Thus by a standard argument we can state there exists a T* = 
T^OIV'oIIl 2 , II Vf/'olU 2 ), such that there exists a unique solution 

V> € C([0,T*);fl' 1 (R 2 )) ni 9 ([0,r);L r (E 2 )) V (q,r) adimssible pair, 

of the Cauchy problem for (4.2). 

Furthermore, again by standard arguments, namely since the energy is conserved for 
every time, we can repeat the same argument for each space-time slab [kT* , (k + 
1)T*}, since T* depends only on the initial datum, and thus we can extend this 
solution globally. Hence theorem 4.2 is proved. □ 

5. COLLISIONAL CASE: THE FRACTIONAL STEP APPROXIMATIONS 

Following the approach in [2] , we want to use the tools presented in the previous 
sections in order to construct a sequence of approximate solutions to system (1.1). 

Definition 5.1. {{p T , J T )} is a sequence of approximate solutions to the system 
(1.1), with locally integrable initial data (p , J ), if (y/fF, A T ) <= Lf oc ([0, T); iJ/ oc (M 2 )) x 
LlMn L toc( R2 ))> P T ■= (vW J T ■= yfW^, satisfy 

• for any test function n G Cg°([0,T) x R 2 ) 

(5.1) f [ p T d t r] + J T -Vr]dxdt+ [ p o v(0)dx = o(l) 

JO JR 2 JR 2 

as r — > 0, 

• for any test function C e Cg°([0, T) x R 2 ; R 2 ) 

(5.2) f [ J T ■d t ( + A T <S)A T : \7( + P(p T ) div (- p T \7V T ■(- J T ■( 
Jo Jm 2 

+ h 2 Vyj r p~?® Vy/y-. VC- -r /o T AdivCdxdt + / Jo ■ C(0)dx = o(l) 

4 Jr2 

as r — > 0; 

• (generalized irrotationality condition) in the sense of distribution 

VA J T = 2VVp^AA T . 
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Our fractional step method is based on an operator splitting method which 
composes the solution operator of the non-collisional problem with the approximate 
solution operator (semi-discrete explicit Euler scheme) of the collisional problem. 
Given a (small) parameter r > 0, we solve a non-collisional QHD problem as in 
Section 4, hence we solve the collisional problem without QHD. At this point in 
order to restart with the non-collisional QHD problem, the main difficulty which 
needs to be solved regards how to update the initial data at each time step. Indeed, 
as remarked in the previous section we are able to solve the non-collisional QHD 
only in the case of Cauchy data compatible with the Schrodinger picture. This 
fact imposes the necessity to reconstruct a wave function at each time step, which 
is made possible by the polar factorization of the wave function computed in the 
previous time step. 

As previously remarked, this method requires special type of initial data (po, Jo), 
namely we assume there exists ipo E H 1 (M 3 ), such that the hydrodynamic initial 
data are given via the Madelung transforms 

Po = IV'ol 2 , Jo = ftlm(V'oVV'o). 

This assumption is physically relevant since it implies the compatibility of our 
solutions to the QHD problem with the wave mechanics approach. The iteration 
procedure can be defined in this following way. First of all, we take r > 0, which 
will be the time mesh unit; therefore we define the approximate solutions in each 
time interval [kr, (k + I)r), for any integer k > 0. 

At the first step, k — 0, we solve the Cauchy problem for the Schrodinger-Poisson 
system 

( ihd t ip T + ^-AiP T = l^l^- 1 ^ + V T tp T 

(5.3) I -AV T = \ip r \ 2 
{ ip T (0) = V'o 

by looking for the restriction in [0, r) of the unique strong solution ip <G C(K; H 1 (M 2 )) 
(see Section 4). Let us define p T := \ip T \ 2 ,J T := ftIm(^ T V7/> T ), then, as shown in 
the previous section, (p T , J T ) is a weak solution to the non-collisional QHD system. 
Let us assume that we know ip T in the space-time slab [(k — l)r, kr) xR 2 , we want 
to set up a recursive method, hence we have to show how to define ip T , p T , J T in 
the strip [kr, (k + 1)t). 

In order to take into account the presence of the collisional term / = J we update 
ip T in t = kr, namely we define ip T (kT+), in such a way that its hydrodynamical 
quantities satisfy 

(5.4) P T (kT+) = p T {kT-), J T {kT+) = (I — T)J T (kT-), 

which means, the hydrodynamical quantities approximately solve (as a first order 
approximation) the collisional problem without QHD. This could be done by means 
of the polar decomposition, described in Section 3. Indeed, denote by 

where clearly (\> T k is a polar factor for the wave function ip T (kT—). Hence a way of 
defining Tp T (kr+) might be the following one 

(5.5) V(kr+) = (^) (1 - r) ^(fcr-), 

in this way the equations (5.4) would be exactly fulfilled, but unfortunately, as we 
will see in the next Section, this way of updating the wave function ip T presents 
some problems, due to the lack of differentiability of the polar factor <p, when we 
try to get a priori bounds for the sequence of approximate solutions. In order to 
circumvent this difficulty we will use approximations for the updating in (5.5) which 
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W> T (fcr+) = W r (/cr-) - i-(t>*k T {kT~) + r k ,r, 



are stable in the energy norm, namely we will apply Lemma A.l in the Appendix, 
with tj) = ip T (kT— ), e = r2 _fe ||?/;o||_H" 1 (R 2 )- In this way we can define 

(5.6) iP T (kT+)=4>, 

by using the wave function tp defined in the Lemma A.l. Therefore we have 

(5.7) p T (kr+) = p T (kT-) 

(5.8) A T (fcr+) = (l-r)A T (fcr-)+i? fc , 

where ||i? fe || L 2 (R 2 ) < T2- fe ||^ ||ffi(R2) and 
(5.9) 

with \\4>*\\ L °° < 1 and 

||r fe , T || L2 < C{T\\Vr{kT-)\\+T2- k \\M\HHm) £ tE$. 

Let us remark that (5.9) will play a fundamental role to prove Strichartz-type es- 
timates for the approximating sequence. 

We then solve again the Schrodinger-Poisson system with initial data ip(Q) = 
ip T (kT+). We define ip T in the time strip [fcr, (k + l)r) as the restriction of the 
Schrodinger-Poisson solution just found in [0,r), furthermore, we define p T := 
\ip T \ 2 , J T := hlm(ip T S7ip T ) as the solution of the non-collisional QHD (4.1). 
With this procedure we can go on every time strip and then construct an approxi- 
mate solution (p T , J T , V T ) of the QHD system. 

Theorem 5.2 (Consistency of the approximate solutions). Let us consider the 
sequence of approximate solutions {{p T ', J T )} T >o constructed via the fractional step 
method. Assume there exists (y/p, A) G L 2 OC ([0, T); iJ^M 2 )) xL 2 oc ([0, T); J 2 oc (R 2 )), 
such that 



(5.10) 
(5.11) 



A T — > A 



m L 2 ([0,T);Hl oc (M. 2 )) 
m L 2 ([0,T);J 2 OC (M 2 )). 



Then the limit function (p, J), where as before J = ^/pA, is a FEWS to the QHD 
system, with Cauchy data (po, Jo)- 

Proof. Let us plug the approximate solutions (p T , J T ) in the weak formulation and 
let C € Cg°([0, T) x K 2 ), then wc have 

/ / J T -dt( + A T <E)A T : V( + P(p T ) div (- p T VV T J T • C 

JO Jr 2 

h 2 f 

+fi 2 VV7 ® : VC - —p T A div (dxdt + J • C(0)dx 

00 r(k+V)T i- 

V / / J T ■ dtC + A T <£) A T : VC + P{p T ) divC - p T W T ■ ( - J T ■ ( 

u-n JkT Jr 2 

h 2 f 
+fi 2 VV7 ® Vy/Jf 7 : VC - —P T A div Cdccdt + / J • C(0)dx 

4 JR2 

r r(k+V)T r- r- 
V / / -J T -Cd2xh+ / J T ((k+l)T-)-(((k+l)T) 

k=Q yJkr Jr 2 Jr 2 



k=0 



-J T {kT)-C(k T )dx 



R 2 



J • C(0)dx 
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= V / / - J T • Qdxdt + V / (J T (fcr-)- J r (fcr))-C(fcr)da 

= V/ / -J T -(dxdt + y^T J T (kr-) ■ ((kr)dx 

k=0 Jkr Jm? k=1 Jm? 

00 „ 

+ J2 / (l-r)V7Mi?fc-C(fcr)dx 

fci^ 2 

00 /-(fe+l)r 
fe=0 ifer ./*» 

00 „ 

+ (1 -t)Y, VF(kr)R k ■ C(fcr)dx 



fc=0' 

o(l)+0(r) 



as r — > 0. 



□ 



6. A PRIORI ESTIMATES AND CONVERGENCE 

In this section we obtain various a priori estimates needed to show the compact- 
ness of the sequence of approximate solutions [p T , J T ) in some appropriate function 
spaces. As we stated in the Theorem 5.2, we wish to prove the strong convergence 
of {^} in Ll c ([0,nHlM 2 )) and {A-} in Lf oc ([0,T); Lf oc (M. 2 )). To achieve 
this goal we use a compactness result in the class of the Aubin-Lions's type lemma, 
due to Rakotoson-Temam [24] (see the Section 2). The plan of this section is the 
following, first of all we get a discrete version of the (dissipative) energy inequality 
for the system (1.1), then we use the Strichartz estimates for Vtp T by means of the 
formula (6.4) below. Consequently via the Strichartz estimates and by using the 
local smoothing results of Theorems 2.8, 2.9, we deduce some further regularity 
properties on the sequence {Vtp T }. This fact will be stated in the next Proposition 
6.7. In this way it is possible to get the regularity properties of the sequence {Vip T } 
which are needed to apply the Theorem 2.11 and hence to get the convergence for 
{^} and {A T }. 

Concerning the energy inequality we notice that for a sufficiently regular solution 
to the QHD system one has 



(6.1) E{t) = - I I |A| 2 dxdt' + £ , 

Jo Jm 2 

where the energy is defined as in (1.7). Now we would like to find a discrete version 
of the energy dissipation (6.1) for the approximate solutions to (1.1). 

Proposition 6.1 (Discrete energy inequality). Let (p T ,J T ) be an approximate so- 
lution to the QHD system, with < r < 1. Then, for all t G [At, (A + l)r) it 
follows 



(6.2) 



N 

E T (t) ||A(fcr-)|| L2(R2) + (1 + T )Eo. 

fe=i 
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Proof. For all k > 1, we have 

E T (kr+) - E T (kr-) = J i|A T (fcr+)| 2 - 1 \A T (kr-)\ 2 dx 

"(-2r + r 2 )|A r (fcr-)| 2 



+ 2(1 - r)A T (fcr-) • R k + \R k \ 2 dx 
<I J (-2t + r 2 )|A T (fcr-)| 2 dx + (1 - r)a|A T (fcT-)| : 

+ ^|i2 fc | 2 + |ifc| 2 da; 
= i J (-2t + r 2 + a - aT)|A r (fcr-)| 2 



Here i?/c denotes the error term as in espression (5.8). If we choose a = t, it follows 

(6.3) ^(fcr+) - £7(fcr-) < -pA T (fcr-)||£ a + | 2 

<_l||A- ( fcr-)|| 2 2 +r2- fe - 1 ||^||H 1 . 

The inequality (6.2) follows by summing up all the terms in (6.3) and by the energy 
conservation in each time strip [fcr, (fc + l)r). □ 

Unfortunately the energy estimates are not sufficient to get the necessary com- 
pactness to show the convergence of the sequence of approximate solutions. Indeed 
from the discrete energy inequality, we only get the weak convergence of V^ 1 " in 
L°°([0, oo); H 1 (R 2 )) 7 and therefore the quadratic terms in (5.2) could exhibit some 
concentrations phenomena in the limit. 

More precisely, from energy inequality we get the sequence {tp T } is uniformly 
bounded in L°°([0, oo); if^R 2 )), hence there exists V € L°°([0, oo); if^R 2 )), such 
that 

V> T ^V inL oo ([0,oo);iJ 1 (R 2 )). 

VF^VP in J L oo ([0,oo);i/ 1 (M 2 )) 
A r ^A inL°°([0,oo);L 2 (R 2 )) 
J^-VPA mL°°([0,oo);L} oc (R 2 )). 

The need to pass into the limit the quadratic expressions leads us to look for 
a priori estimates in stronger norms. The link with Schrodingcr equation brings 
naturally into this search the Strichartz-type estimates and the following results are 
concerned with these estimates. However they are not an immediate consequence 
of the Strichartz estimates for the Schrodinger equation since the wave function 
obtained in the updating procedure implement in the fractional step in general 
does not solve any Schodinger equation. 

The following Remark 6.3 will clarify in detail why we used the result in Lemma 
A.l to update the wave function in place of using the exact factorization property 
as in equation (5.5). 



Therefore we get 
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Lemma 6.2. Let ip T be the wave function defined by the fractional step method 
(see Section 5), and let t G [Nt, (N + 1)t). Then we have 

N 

(6.4) VV T (t) - U(t)Vi> -i^Yl W - fcT ) (^A^(fcr-)) 



fc=i 



i / U(t- s)F{s)ds + y^Uit- kr)r T k , 
Jo fc=i 



w/iere as before U{t) is the free Schrodinger group, 

(6-5) H$fclU°°(R3) < i, lkfclU 2 (K 2 ) < ^II^oIIhmr 2 ) 

and F = V^Y'V + V T V)- 

Remark 6.3. Now we can see more in details why the updating step in (5.6) has 
been defined through the Lemma A.l. Indeed, one could possibly try to use the 
exact polar factorization as in the Lemma 3.1, namely ip T (kT+) := 0( 1_T ) ^fp, where 
4>, y/p are respectively the unitary factor and the amplitude of ip T (kT— ) (note that 
</) is uniquely determined ^ypdx— a.e., and \<p\ — 1, y^pdx-a.e.), as we anticipated in 
Section 5, see equation (5.5). In this case for the hydrodynamic quantities we get 
the exact updating formulas 

W(kr+)\ 2 = \V{kr)\ 2 

hlm(i^Wtp T )(kT+) = (1 - T)ftIm(V^VV> T )(fcT-). 

In order to obtain Strichartz-type estimates for the approximating sequence we need 
to compute VV> T , the gradient of the approximated wave function. This calculation 
involves at time t, via the Duhamel's formula, the updates done at each kr < t: 

VV> T (t) =U{t - Nt)<t t Nt U(t) . . . a T T U(r)Vi;o 

_ iLy(t - Nt)<j t Nt A t (Nt-) + ...- i T jU(t -Nt)... U{t)4> t t A t (r-) 

i I U(t- s)F(s)ds - iU(t - Nt)g t Nt / U(Nt - s)F(s)ds 

J Nt J(N-1)t 

+ ...-iU(t- Nt)o- t Nt U{t)...<t t t f U{t - s)F(s)ds, 

Jo 

where aJ, T = ((f> k ~ T )~ T , <p kT being the polar factor of %p T (kT~). Now, to use the 
group properties of the free Schrodinger operator, we need to pack all the U(-Ys 
and all the ct t 's. The estimates of these commutators [[/(•), cr T ] are in general not 
true for a T kT which are only in L°°(R 2 ) (sec [27]). 

Proof. Since ip T is solution of the Schrodinger-Poisson system in the space-time 
slab [Nt, (JV + 1)t) x I 2 , then we can write 

(6.6) W^ T (t) = U(t- NT)Vip T {NT+)-i f U(t- s)F(s)ds, 

JNt 

where F is defined in the statement of the Lemma 6.2. Now there exists a piecewise 
smooth function 9n, as specified in the proof of the Lemma A.l, such that 

and furthermore all the estimates of the Lemma A.l hold, with ip — tp T (NT—), 
ip = ■0 r (A r r+) and e = 2 _Ar r||')/'o||ff 1 (R 2 )- Therefore, we have 

(6.7) Vip r (Nt+) = V^ t {Nt~) - i^ 1 - T)eN A T (NT-) + r T N , 
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where Hr^H^ < tH^o || jt 1 - By plugging (6.7) into (6.6) we deduce 
Vip T (t) = U(t - Nt)ViI> t (Nt-) - i T jU{t - NT){e l{1 - T)6N k T (NT-)) 



+ U(t - Nr)r T N -if U(t- s)F(s)ds. 

JNt 



Let us iterate this formula, repeating the same procedure for Vi(j t (Nt—), then (6.4) 
holds. □ 

At this point we can use the formula (6.4) to obtain Strichartz estimates for 
VV> T , by applying the Theorem 2.6 to each of the terms in (6.4). After some 
computations, the following result holds. 

Proposition 6.4 (Strichartz estimates for Vip T ). For all T > 0, let ip T be the 

approximating wave function, then one has 

(6-8) l|V<A T || ifi . ([0 ,T]x R 2) <C(eI,\\ Po \\ lhwl2) ,T) 

for each admissible pair of exponents (q,r). 

Proof. First of all, let us prove the Proposition 6.4, for a small time < T\ < T 
and let (q, r) be an admissible pair of exponents. We will choose T\ > later. Let 
N be a positive integer such that Ti < Nr. By applying the Theorem 2.6 to the 
formula (6.4), we get 

L«LJ([0,Ti]xR 2 ) 

JV 

\\u(t - kr) (e t{1 - T)9k K T {kT-) 

k=i 

N 
k=l 

f U{t-s)F(s)ds 
Jo 



T 

h 



£?£E([0,Ti]xR 2 ) 



L'LJ ([0,71] xR 2 ) 



=: A + B + C + D. 



Now we estimate term by term the above expression. 
The estimate of A is straightforward, since 

||^(i)W>o||L t 9 Lj([0,T 1 ]xR2) iS ||V-(/>o||l2(R2)- 

The estimate of B follows from 



N 

(6.9) ^ Yl - kT ) (e i(1_T)flfc A T (fcr-)) 



fc=i 



L?Lj([0,Ti]xR 2 ) 
N 

fe=i 



The term C can be estimated in a similar way, namely 



N 



N 



(6.io) J2 IW* - ^yihtLr < E MWw) Z TiUo\\H H m- 



k=l 



k=l 



The last term is a little bit more tricky to estimate. First of all we decompose F 
into three terms, F = F ± + F 2 + F 3 , where F 1 = V(|^ r | p_1 ^ r ), F 2 = VV T Tp T and 
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F3 = V T Vip T . By the Strichartz estimates (Theorem 2.6), we have 



(6.11) 



/ U(t-s)F(s)di 
Jo 





<ll^i 



Z^LjQO/r^xR 2 ) 



l*3| 



where (qi,Vi) are pairs of admissible exponents. Now, we summarize the previous 
estimates by using (6.4) in the following way 

(6.12) ||VVlso ([0iTl]xK2) < ||VVo|U 2 ( r 2 ) +T,eI +Tr||V^|| P SO([0 , Tl]xR2) 

S°([0,Ti]xR 2 ) 



< (1 +T)E* +T? +TMMh m m T \\so ao , Tl] ^y 



Lemma 6.5. There exist T^Eq, \\^o\\ l ^),T) > andd(E , \\i>o\\ L '(R*),T) > 0, 
independent on t, such that 

(6-13) IIW> T |Uo ([0i f]x R2) < CtiEo, \\M\W),T) 

for allO <f<T x (£ , ||V>o||.l 2 (r 2 ))- 

Let us recall that E T (t) = E n and ||?/; T (i)|| i 2( R 2) = ||V'o||l 2 (r 2 ), hence we are in 
the situation in which we can repeat our argument on every time interval of length 
Ti, always depending on the same parameters Eo, \\iPo\\l 2 - Hence this yields to the 
following inequality on [0,T] 



3.14) \\vrw S°([0,T]xR 2 ) 



<C\(E q ,\\Mli,T) 



T 



+ 1) =C(\\Ml*,E ,T). 



□ 



Proof of the Lemma 6.5. Let us consider the non-trivial case \\ipo || L i > 0. Assume 
that X € (0, 00) satisfies 

(6.15) X < A + f iX + XX p = cf>(X), 

with p > 1, A > and for all < p < 1, A > 0. Let X« be such that <j>'(X*) = 1, 
1 

namely X* = {^^j" \ hence one has <f>(X*) < X* each time the following 

' \ 1 \ (1 - p)^ 



inequality is satisfied 
(6.16) 



pi 



> A. 



Therefore the convexity of <f> implies that, if the condition (6.16) holds, there exist 
two roots X±, X + (p, X,A) > X* > X_(fi,X,A), to the equation <j>{X) = X. It 

then follows either < X < X_, or X > X + . In our case p = 7\ 1/2 || Mollis, A = Tf , 

1 /2 

A=(l+T)Eo', hence we assume 

1 



v = Tl /2 \\Mh< 



2 

— P 
P P-1 



[2— (1 + T)E* 



1/2 



p-i 
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Therefore we choose 



(6.17) 



Ti := min 



(2||V'oIIl 2 )- 2 , 



2^t(1 + T)£ 1/2 _ 



Clearly we cannot have 

"l-Woll^ 



z* = 



< X + < ||VY> r || i go([o j T 1 ]xIR 2 )' 



since we get a contradiction as T\ — >• 0, hence 



(6.18) 



||W T 



IS°([0,Ti]xR 2 ) 



□ 



Corollary 6.6. For all T > 0, let \f(f , A T be the approximating sequence defined 
by the fractional step method, then 

i 

(6.19) \\^VW\\LlL-([0,T]xm) + I|A T ||l?L-([0,T]xR2) < C(E^,\\p 1| L i (R 2) , T), 

for each admissible pair of exponents (q,r). 

Unfortunately this is not enough to achieve the convergence of the quadratic 
terms. We need some additional compactness estimates on the sequence {Vtp T } in 
order to apply Theorem 2.11. In particular we need some tightness and regularity 
properties on the sequence {Vf/> T }, therefore we apply some results concerning local 
smoothing due to Vega [29] and Constantin, Saut [9]. 

Proposition 6.7 (Local smoothing for Vtp T ). For all T > 0, let ip T be defined as 
in the previous section, then it follows 

(6-20) IIW1l W ];<?(R 2 )) " C(^o,T,||p |UO- 

Proof. Using the Strichartz estimates obtained above, we can apply the Theorems 

2.8, 2.9 about local smoothing. Indeed, by using again the formula (6.4) it follows 



II V ^ T Hl 2 ([0,T]:</ c 2 (R2)) ^I|V^o||l 2 (R 2 ) 

N 

+ T Y J W(kr-)\\ L ^ ) 
%=i 

N 

+ T^||V?/v i J| i 2 (R 2 ) 



L 2 (R 2 ) 



k=l 

N 

+ E || V ^ " WfcT-) + Uk nk - A T (kr-)) 

+ l|-F||Li([0,T];L 2 (R 2 ))- 

□ 

Furthermore, since again we have the formula (6.4); for t € [Nt, (N + l)r), 
(6.21) V?A T (t) = U(t)V4> Q -i [ U(t-s)F(s)ds 

N N 

'h 







- 4 E u $ - kr ) [<« AT (^-)] + E u (* kT >^ 



k=l k=l 
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holds (see Chapter 6), we can exploit the Strichartz estimates to obtain a bound 
like the one in (4.10), 

(6-22) ||VV T |U«i S ([o,r|xR a ) <C(T,\\ip \\ L 2,E ), 

for every admissible pair (q,r) of Schrodinger exponents in R 2 . Hence, in order 
to obtain the sufficient compactness properties on the sequence of approximate 
solutions {V-0 T }, it only remains to recover the smoothing estimates for Vip T . But 
this is straightforward since by theorems 2.8, 2.9 we have 

(6.23) 

L*([0,T\;Hl£(m?)) ~ IIW'o||l,2 + |||V>T 1 ^ ^ II L\Ll + 11 VVll>\\ L i L 2 +\\ Wlp\\ L i L 2 

< HVVolU- + ^||V||^ i? ( P - 1 )||vviL rL? 

+ ^11^11^11^11 & +r||W|| L oo iS ||V|| ^ 

< \\vML'+T^\m% Li \m\^~ i 

+ r^||V-||z.^z. £ ||V^||^^^^ + ^11 VV-ll^^^gll-011^^.11 V-0llf^ iS , 

where (<?i, ri) is a Schodinger admissible pair in K 2 , r 2 is such that 2 < r 2 (p — 1) < 
oo, r is such that 1^(0) <G U and 2 < p < oo. Note that in the second inequality we 
used Holder's inequality, while in the third one the Gagliardo-Nirenberg-Sobolev 
inequality. Hence, by the estimate (6.22) we can say 

(6.24) IIW1 LW];<?(R2)) < C(T,Uo\\ L2 ,E ). 



Hence, we can apply the theorem 2.11 by Rakotoson, Temam; this means there 
exists ip G L' x> ([0,T];i2" 1 (R 2 )) n £ 2 ( [0, T] ; iff/j 2 (M 2 ) ) such that 

(6.25) V T ^-0 inl^CT];^ 1 ^ 2 )) 

strongly. Finally, by applying the consistency theorem 5.2 we can state that (p, J) 
defined by (p, J) := (\tp\ 2 , hlm(ipVip)) is a finite energy weak solution to (1.1), (1-2). 

Appendix A. Polar Decomposition 

In this Section we will get into details of the main properties of polar factoriza- 
tion of a wave function, in particular we will prove Lemma 3.1. 

Proof of Lemma 3.1. Let us consider a sequence {ip n } C C°°(R 2 ), ij) n — > -0 in 
_ff 1 (R 2 ), define the unitary factor 

<"> 

Then, there exists <p <= L°°(R 2 ) such that <f> in L°°(R 2 ) and V^„ -> in 

L 2 (R 2 ), hence 

Re (4^Vil> n ) Re (0Vr/>) in L 2 (R 2 ). 
Since by (A.l), one has 

Re(^VV>„) = V|V„| 

a.e. in R 2 , it follows 

VV^^Re(^W) ini 2 (R 2 ). 
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Moreover, one has Vy/p^" — 1 V^/p in L 2 (R 2 ), therefore 

VVP = Re (0Vr/>) , 

where ^ is a unitary factor of ^i>. 

The identity (3.4) follows immediately from 

(A.2) fi 2 Rc(a^9fcV') = fi 2 Rc((<?!)a^)(09feV)) = h 2 Re((j>d 3 ^)Re($d k rb) 

- ft 2 Im(0a^)Im(^a fe V) = h 2 d is /pd ks fp + A^A^. 
This identity is valid for any ip £ H 1 (M 3 ) since Vip = a.e. in the nodal region 

{V = o}. 

Now we prove (3.6). Let us take a sequence ip n — > ip strongly in H 1 (R 2 ), and con- 
sider Vy/fa — Re (0nV^ n ), A„ := ftlm (0 n VV'n)- As before, <j> n in L°°(R 2 ), 
with (/> a polar factor for t/>; then V^/p^ — 1 V^p, Re (^Vi/'n) — Re ((f>Vip), and 
V^/p = Re(^VV'). Moreover, A„ := Mm(^Vf„) film (<£v^>) =: A. To up- 
grade the weak convergence into the strong one, simply notice that by (3.4) one 
has 

fi 2 ||W||| 2 = ft 2 ||VVp||| 2 + ||A||| 2 < hminf (fi 2 || V^H^ + ||A„|| 2 2 ) 

n— >oo 

= liminf ft 2 ||V^||| 2 =H 2 \\VtP\\ 2 l2 . 

n—>oo 

Finally to prove identity (3.5) it suffices to notice that we can write 

nv^AVtp = h{4>vi>) a (0VVO, 

where <j> € L°°(R 2 ) is the polar factor of ip such that V-^/p = Re((pVip), A := 
fiIm(</>V^>), as in the Lemma 3.1. By splitting (pVip into its real and imaginary 
part, we get the identity (3.5). □ 



Now we state a technical lemma which will be used in the next sections. It summa- 
rizes the results of this section we use later on and will be handy for the application 
to the fractional step method. 

Lemma A.l. Let ip G H 1 (R 2 ), and let r, e > be two arbitrary (small) real 
numbers. Then there exists ip g H 1 (M. 2 ) such that 

P = P 

A = (l-i-)A + r e , 

where ^fp := \ip\,y/~p :— \tp\,A :— hlm(tpVip), A := hlm((pVip), (p,<p are polar factors 
for ip, ip, respectively, and 

|ke||L2( K 2) < e. 

Furthermore we have 

(A.3) Vip = Vip-i^(p*A + r £iT , 

w/iere ||</>*||l°°(r 2 ) < 1 and 

(A.4) llr^ll^^) < C(r||V^|U2(R2) + £)■ 

Proof. Consider a sequence {ip n } C C°°(R 2 ) converging to ip m H 1 (R 2 ), and define 

, . . f TT^TT if^AnW 7^0 
if 1p n (x) = 
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as polar factors for the wave functions ip n . Since ip n € C°°, then <p n is piecewise 
smooth, and 0„ := {i £ I 2 : |^ n (a;)| > 0} is an open set, with smooth boundary. 
Therefore we can say there exists a function n : fl n — > [0, 2ir), piecewise smooth in 
fl n and 

4> n (x) — e l6n{x \ for any i e 0„. 

Moreover, by the previous Lemma, we have <f> n <p in i°°(M 2 ), where 4> is a polar 
factor of V, and A„ := SIm(^) n Vi/) n ) -> A := Mm(^V^) in L 2 (]R 2 ). Thus there 
exists n€N such that 

\\lp ~ VVl||/fi(IR2) + ||A - A„|| L 2( R 2) < E. 

Now we can define 

(A.5) i,:=e^-^^. 

Let us remark that (A.5) is a good definition for tp, since outside fi„, where 6 n is 
not defined, we have y/pn~ = 0, and hence also ip n — 0. Furthermore, 

n 



+ (Vi - W) - 2^e i ( 1 -) e »(A„ - A) 
n 

and thus we obtain (A. 3) and (A. 4), with r £jT given by 

£ H^' 1 ) W„ + (Vi - V^) - i^-^(A n A) 




(A.6) A = Hmle-'C-^-V^) - (1 - r)A„ = (1 - r)A + (1 - r)(A„ - A) 
and clearly r e := (1 — r)(A„ — A) has L 2 -norm less than e by assumption. □ 

Acknowledgement. We wish to thank prof. R. Danchin for having drawn our 
attention to an important technical remark needed to ensure the validity of the 
identity (A.2). 
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